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Abstract. The collective pairing hamiltonian is obtained in the framework of the generator coordinate 
method in the gaussian overlap approximation with a slightly modified BCS function used as a generator 
function. The collective variable a, measuring the monopole moment of the pairing field, and the gauge 
transformation angle <f> are chosen as generator coordinates. The vibrational ground states are calculated 
by diagonalisation of the collective pairing hamiltonian in the harmonic oscillator basis. 

PACS. 21.30.-x Nuclear forces - 21.60. Ev Collective models 



1 Introduction 

The collective pairing hamiltonian was originally intro- 
duced in Ref. [1], where the intrinsic deformation of the 
pairing field a, related to the BCS gap parameter A, and 
the gauge transformation angle (j> were used as collective 
variables. The pairing hamiltonian was then derived in 
framework of the cranking approximation. A more gen- 
eral method which determines wave functions and the col- 
lective hamiltonian is the Generator Coordinate Method 
(GCM). The GCM derivation of the collective hamiltonian 
with the BCS functions as the generator functions were al- 
ready performed in Ref. [2-4], where the Gaussian Over- 
lap Approximation (GOA) of the generator wave functions 
was used. The collective pairing Schrodinger equation was 
constructed using the single-particle plus monopole pair- 
ing hamiltonian. 

The aim of this paper is to derive collective vibrations 
of the pairing field in GCM+GOA approach for the hamil- 
tonian based on the single-particle Nilsson potential [5] 
and the ^-pairing interaction [6-8] . The latter requires in- 
troducing a suitable collective variable different from the 
monopole pairing gap used in [2]. Here we have used the 
same generator coordinate a = u kVk as in Ref. [1] , 
which as suggested in Ref. [9] is a natural choice for this 
type of collective motion. 

The paper is organized as follows: In Sect. 2 we dis- 
cuss the properties of the 5-pairing interaction. In Sect. 
3 the collective coordinates and the form of the collective 
pairing hamiltonian are introduced. The formulae for the 
GCM+GOA mass parameters are compared to the crank- 
ing ones. The method of diagonalising the collective pair- 
ing hamiltonian is described in Sect. 4. Sect. 5 contains 
numerical results and Sect. 6 conclusions. In Appendix 
A we discuss the multipole expansion of the <5-force and 



demonstrate validity of our choice of a as the collective 
coordinate. 



2 5-pairing forces 

The nuclear mean-field hamiltonian with the residual pair- 
ing interaction can be written as 



H = H. 



S . J J . 



H 



pair i 



where 



H s . p .=J2( k \Mk)(clc k + clc h ), 

Vkkll c k C k C l Cl ■ 



k>0 



(1) 
(2) 
(3) 



k,l>0 



The summation runs over the eigenstates \k) of the single- 
particle hamiltonian h. The antisymmetrized matrix ele- 
ment of the pairing interaction V^u m the pairing hamil- 
tonian (3) is given by the following expression 



Vkkll 



cr 1 a 2 

V T (r 1 ,er 1 ;r 2 ,(T2) 
[$i(ri,<Ti)$i(r 2 ,o-2) - $i(rio-i)$i(r 2 ,(T 2 )] . (4) 



Here <Pk(r,cr) is the single-particle eigenfunction of h in 
the space (r) and spin (er) representation and <Pj,(r,cr) is 
its time reversal counterpart. V r (rx, <T\\ r 2 , (T%) is the im- 
pairing force [6] 
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V T (rx, <Tx;r 2 , er 2 ) = VJ- — — S(n - r 2 ) 



(5) 
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In the following we consider the proton-proton (p-p) and 
neutron-neutron (n-n) part of the interaction only. Taking 
this into account one obtains: 



where 



Gki = V kM = V / d 6 rp k {r) P i{r) , 



Pk(r) = \$ k {r)Y 



(6) 



(7) 



and Vq is the pairing interaction strength as adjusted in 
Ref. [10]. The pairing gap equations become 



4 = ^ Gki 



A, 



l>0 



(8) 



where e; is the single particle energy and the Fermi level 
A is determined from the particle number equation 



efe - A 



k>o \ \J (e k - A) 2 + A k 



(9) 



The coupled set of Eqs. (8) and (9) is solved numerically 
by the appropriate iteration procedure. 



3 Collective pairing hamiltonian 

Following Ref. [11] we construct the monopole pairing op- 
erator (see Appendix) which characterizes the amplitude 
of excitations connected with vibrations of the pair den- 
sity 

A=W{e-^clc\ + e^c kCk ) . (10) 

fc>0 

The pair condensate can be described by the mean value 
of the operator (10) in the BCS like state 



fc>0 



where N is the number of particles, cf) is the gauge angle 
and the average pairing gap is given by the expectation 
value of the operator (10) 



a = (a4>\A\a4>) = ^ u k v k . 



(12) 



fe>0 



Using (11) as a generator function and a, as generator 
coordinates and following the steps of Ref. [2] the collec- 
tive hamiltonian 

K 2 d n , d 

Hon = - —y/Aet laa M aa — 

2vdet7 QCe oa oa 



--% M 



^ d<j> 2 



Im(a</>| ■£rH\a<f>) Q 
—ih tt- 



+V(a), 



(13) 



is derived. In Eq. (13) 7 aa , are related to the widths of 
the Gaussian overlap and -M~^, ■M'll are the components 
of the inverse mass tensor. The quantities appearing in Eq. 
(13) are analogous to the general expressions obtained for 
the A and 4> coordinates in the monopole pairing case [2] . 
V(a) is the collective pairing potential equal to 



V{a) = {act>\H\a4>) - £ , 



(14) 



where £0 is the so-called zero-point energy. The collective 
pairing hamiltonian (13) is hermitian with the Jacobian 



dr = ^/7q Q 700 ada < 



(15) 



It means that the eigenfunctions of (13) should be orthog- 
onal with the above measure. 

The mean value of the hamiltonian (1) is evaluated 
using the constraints N = const and a = const which 
leads to minimization of the average value of the operator 

H' = H — X(N — (ac/)\N\a<j>)) - £(A - (a(j>\A\a<f>)) , (16) 

calculated in the BCS state (11). Here A and £ are La- 
grange multipliers. 

The expectation value of the BCS hamiltonian in Eq. 
(14) is given by 

Ebcs = (a<j>\H\a(f>) = 

= 2E e ^-^E#-E<W< 



fc>0 



2 ^ E k 

fe>0 fc>0 



where E k — y/ (e k — A) 2 + A\ is the quasiparticle energy. 

The final expressions for the nonvanishing components 
of the metric and mass tensors are 



7aa = (j2 E 

k 



(e k - A) 2 
16E£ 



A 2 

E^ k 

k k 

k 



laa ) 



where a is equal to 



E 



(e k - A) s 



(18) 

(19) 
(20) 

(21) 
(22) 



The zero-point energy appearing in Eq. (14) reads 



(e fc - A)' 
16^ 3 



/ 



.(23) 
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The first term in Eq. (23) represents the collective cor- 
relations in the ground state energy whereas the second 
corresponds to the approximate particle number projec- 
tion [4]. 

An alternative method of treating collective vibrations 
is the cranking model [12, 13]. The mass parameters in this 
case read 

k 



J l aa 



16EI 



A 2 



(25) 



and the collective potential is usually assumed to be equal 
to E BCS (17). 



4 Diagonalisation of the collective pairing 
hamiltonian 

Having found the collective potential, the mass and the 
metric tensors, we can evaluate the collective hamilto- 
nian numerically. Only the vibrational spectra will be con- 
structed here as the quasirotational states correspond to 
the bands built from the ground states function of the 
neighbouring even- even isotopes (or isotones). In the first 
step we perform a transformation from the a coordinate 
to a new variable x in which the mass parameter M xx is 
nearly constant 



ox 



const . 



(26) 



The mass parameter M aa is a rapidly decreasing function 
of a and can be approximated by the function 



M aa = 



(a + ao) 2 



(27) 



where ao and b parameters are chosen to approximate the 
mass parameter in the best way. Variables x and a are 
connected through the equation 




2 4 6 8 10 12 14 16 

a 

Fig. 1. Pairing potentials: the BCS energy (dashed line), the 
GCM collective potential (solid line) and the particle number 
corrected BCS energy Ebcs — £q (thin solid line) as functions 
of the collective coordinate a for protons in 148 Ce. Zero value 
on the abscissa corresponds to the minimum of the BCS energy. 



These eigenstates which satisfy the usual normalization 
conditions 



ipitpjdx = dij 



(31) 



are employed to calculate the matrix elements of the col- 
lective hamiltonian (Eq. 13). The digonalization was per- 
formed in a basis of N max = 19 oscillator shells. 



x = ^/b/Jiln ( 1 + 



a 



(28) 



where n is an arbitrary constant. The basis states used 
to diagonalise the hamiltonian (13) are generated by the 
harmonic oscillator hamiltonian 



h d 2 1 , , 22 



with frequency u> determined from the plateau condition 
of the energy dE/doj 0, where E is the ground state 
energy of the collective hamiltonian (13). In our case only 
the even eigenstates are picked (see [2]) 



H B ^t = {i + , 



0,2,4,... 



(30) 



5 Results 

We have considered a few nuclei in the rare earth region. 
The A = 165 parameter set of the single particle Nilsson 
potential is used [14]. The 5-pairing strengths as obtained 
in Ref. [10] are 



^neutrons = 2 30MeVfm 3 , 

protons = 240MeVfm 3 . 



(32) 



The fit of the Vq parameters was done for the pairing 
window with 2\/15 n levels closest to the Fermi energy 
where n = N for neutrons and n = Z for protons [14]. 
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Fig. 2. Components of GCM and the cranking mass parame- 
ters for protons in 148 Ce. The M aa components are multiplied 
by the factor 100. The components M<f,<p are in natural units. 



Figure 1 shows the pairing potentials for BCS and 
GCM models. The dashed line represents the pure BCS 
energy (17) as a function of the collective coordinate a 
whereas the solid line corresponds to the GCM collective 
potential (14). The particle number corrected BCS en- 
ergy, i.e., the one with extracted £q energy is also shown 
(thin solid line). Zero value on the abscissa corresponds 
to the minimum of the BCS energy. The particle number 
projected energy is about 2 MeV deeper than the BCS en- 
ergy minimum and its position is shifted in the direction 
of larger a. Analogously, the GCM ground state energy 
which additionally contains the zero point energy corre- 
sponding to the pure a vibrations has a minimum shifted 
by a similar amount towards larger a. 

The M aa and components of the mass tensor 

for both GCM and the cranking model are shown in Fig. 
2. The ratio of the a component of the cranking mass 
(dashed line) to the GCM mass (solid line) is close to 2/3 
in the vicinity of the minimum of the GCM potential. The 
decrease of the collective mass with a (or A) has a sig- 
nificant influence on the spontaneous fission half lives of 
heavy nuclei, as shown for A collective coordinate, as well 
as on the height of the fission barriers (see e.g., [15]). Sim- 
ilarly in case of the collective Bohr model calculations [16, 
17] the inclusion of coupling of quadrupole vibrations im- 
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Fig. 3. The collective potential V(a) (solid line), the GCM 
mass parameter M aa (a) (dotted line) and the probability den- 
sity P = |<?o(a)| 2 ^W700 (thin solid line) of the ground state 
for protons in 148 Ce. The ground state energy Eq is marked. 
The arrows indicate the positions of the minimum of the po- 
tential and the maximal value of the probability distribution 
on the a axis. 



proves significantly the agreement with the experimental 
data. 



In Fig. 3 we show the collective potential V(a) (solid line), 
the GCM mass parameter M aa (a) (dotted line) and the 
probability density P = |^o(a)| 2 ^)W70^ of the ground 
state (thin solid line). The short line segment in the mid- 
dle of the figure marks the position of the ground state 
energy. The arrows indicate the positions of the minimum 
of the potential (bottom scale) and the maximal value of 
the probability distribution (top scale). The a value of 
the equilibrium is 9.3 and the most probable value is 7.6. 
The latter is shifted towards smaller a values that im- 
plies the increase of the M aa by 1.5 factor on average. 
This behaviour is common for the majority of nuclei in 
the considered region. As shown in Ref. [16] in case of 
monopolc pairing the use of the most probable value of 
the gap parameter (instead of the BCS gap) leads to a 
considerable improvement of the predictive power of the 
Bohr hamiltonian. 
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6 Conclusions 

The first excited collective pairing vibrational states for 
even-even nuclei in rare earth region appear at energies 
close to 2.5 MeV for protons and 4.5 MeV for neutrons. 
Usually they have higher energies than two quasiparticlc 
excitations and consequently one has to include these cor- 
relations in the ground state, only. 

Since the pairing vibrations are strongly coupled with 
shape degrees of freedom of the nucleus it is hard to com- 
pare the results (which are similar to that obtained with 
the coordinate A) with experimental data. 

Finally, it should be emphasized that the role of the 
collective coordinate a is analogous to the A coordinate 
considered earlier. 



We are indebted to Jerzy Matyjasek and Paul Stevenson for 
careful reading of the manuscript and valuable comments. 



A Multipole expansion of 5(r)-force 

In the following we consider the short range interaction 
operating between equivalent nucleons (T = 1, protons or 
neutrons) 

V 12 (S) = V(r 1 ,<T 1 ;r 2 ,(T2) 

= V Q 1 - (T ^ eT2 ^(r 1 -r 2 ), (33) 

equal to Vo<5 3 (r) in case of the spin singlet (S = 0) and to 
for the spin triplet (S = 1). 

One can use the following representation of 5(r) 

S 3 (rx - r 2 ) = —^—6(ri - r 2 )5(cos6'i - cos6<2)<5(0i - 4> 2 ) ■ 

nr 2 

(34) 

This form is very useful when separating the angular and 
radial parts. Using the relation 

^(cosfli — cos8 2 )5((j)i — (j) 2 ) = 

= E E Y km (di , MYkm (0 2 , 4> 2 )(35) 

K 71 M 

and angular momentum algebra (see Talmi [18]) one gets 
the following energy shift of the identical (T = 1) particles 
state \jij 2 J), 



V S (J) = V F R ( ni hn 2 l 2 ) I(jij 2 J) . 



(36) 



where 



F R { ni hn 2 l 2 ) = ±-j ±Rl ih ( r )Rl 2l2 (r)dr (37) 



and 



iijihJ) 



(2ji + l)(2j 2 + l) ( Ji j 2 J 



1 + S nin2 6 hl2 \, 1/2 -1/2 0, 

(Zi + h + J even) . (38) 



Now we show that the 5 interaction is separable in the 
particle-particle channel. To do this we define the creation 
operators of a pair of particles coupled to the angular mo- 
mentum JM 

P JM = 7| E (ji^iJ 2 m 2 |JM)4 1 cL 2 (39) 



and its Hermite conjugate 



JM 



P 



1M 



(40) 



Here creates a nucleon in (nljm) single particle state. 
The (5-interaction can be written in second quantized form 



JM ■ 



(41) 



J, even 



M 



Note that Vs(J) does not depend on the third component 
M of the total angular momentum of a pair. 

Since the pairing force is the interaction in the particle- 
particle channel with the third angular momentum com- 
ponent M = 0, the operator P JM defined in Eq. (39) sim- 
plifies to Pj = Pj_m=o an d the summation in (39) runs 
over m = mi — —m 2 only. 

The part of the interaction corresponding to J = is 
the classical pairing or better to say monopole-pairing in- 
teraction. The next component J = 2 forms the quadrupole- 
pairing force and so on. The monopolc pairing term reads 



^(o) E 



t t 

C m C m i c m > 



(42) 



The operator c^ creates a particle in the time reversal 
state (m) which is defined by 

{-iy- m d m . (43) 



c - 



The quadrupole (J = 2) part of Eq. (41) reads 

V i( 2 ) E C? 1 ™^ 2 - 777|20)(jim'j2 - m'ftOytnC^Crn'Cm' , 
m.m' 

(44) 

etc. Each term in Eq. (41) separates in a similar manner. 
The strength Vs(J) of each multipole is different and de- 
pends on Vq and the quantum numbers (nihji), [n 2 l 2 j 2 ) 
and J (see Eq. (36)). 

It is easy to show that for e.g., j\ = j 2 = 9/2 the ratio 
V s (2)/Vs{0) w 0.24, Vs{4)/V s {0) w 0.12 etc.. Therefore, 
the state J = has the largest energy shift towards the 
lowest energies (assuming Vq is a positive constant). 

In analogy to Eq. (10) one can define the following 
operators (only even J values are allowed; see eq. (41)) 



(e- 2 ^Pj + e 2 ^Pj) , 



J = 0,2, 4,... 



(45) 



The operators Aj describe different multipolarities of the 
pairing field. Their average values are the multipole defor- 
mations of the pairing field. The mean pairing field Hamil- 
tonian equivalent to the S force expressed in terms of the 
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field operators Aj reads 

V 12 (6) = Yl V s (J)ajAj, (46) 

J , even 

where 

= (Aj) (47) 

is the multipole deformation of the pairing field. The term 
Aq corresponds to the operator A used in the present pa- 
per (see Eq.(lO)) and describes the monopole type defor- 
mation and the parameter ao corresponds to the collective 
deformation parameter a used in the text 

a ~ a ■ (48) 
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